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SHOCK-TURBULENCE INTERACTION AND THE GENERATION OF NOISE 1 

By H. S. Ribnbr 


SUMMARY 

The interaction of a convected field of turbulence with a 
shock wave has been analyzed to yield the modified turbulence , 
entropy spottiness , and noise generated downstream of the 
shock . This analysis generalizes the results of Technical 
Report 1164, which apply to a single spectrum component, to 
give the sh o ck-inter action effects of a complete turbulence field . 
The previous report solved the basic gas-dynamic problem, 
and the present report has added the necessary spectrum analysis . 

Formulas for spectra and correlations have been obtained 
and numerical calculations have been carried out to yield 
curves of root-mean-square velocity components, temperature, 
pressure , and noise in decibels against Mach number for the 
Mach number range of 1 to °° ; both isotropic and strongly axisym- 
metric {lateral p e rturb a tio ns flo ngitudinal perturbations **S6/1) 
initial turbulence have been treated . It was found that in 
either case initial turbulence with a longitudinal component 
of 01 percent of stream velocity wcruld yield a noise pressure 
level of about 120 decibels; the value of lateral component had 
relatively little effect 

The present results are applicable quardilatively to flow in 
duds or channels containing normal shocks; they are presumed 
to provide a qualitative guide to the generation of noise by the 
shock structure in a supersonic free jet . 

INTRODUCTION 

The propulsion of aircraft by means of jets gives rise to 
intonse noise as an unfortunate byproduct. Programs of 
noise abatement are under wa y, but at present they are 
largely empirical; even with the general guide provided 
by Ligh thill's theory (ref, 1), the understanding of the 
mechanisms of noise generation is far from complete. It 
appears from both experimental and theoretical evidence, 
however, that the interaction of turbulence with shock waves 
must often play a part. On the theoretical side, the genera- 
tion of noise by such interaction is pointed out independently 
in references 2 and 3. The shock-turbulence interaction 
was found to produce, in addition to the noise, an entropy 
“spottiness” aft of the shock (manifested as a temperature 
and density spottiness at constant pressure, ref. 2) . 

Turbulence, entropy spottiness, and noise (pressure 
fluctuations) are examples of the three fundamental modes 
of small disturbance perturbation of a gas (refs. 4 and 5) : 
more specifically, the categories are vorticity mode, entropy 
mode, and sound mode. The vorticity mode (turbulence) 


and the entropy mode are essentially frozen” patterns 
(to use Kov^sznay's term) that are convected by the main 
flow; the sound mode, however, consists of waves that pro- 
pagate in various directions in addition to being convected. 

To the first order in the perturbation velocity, there is no 
tendency for the modes to interact or for an isolated mode 
to spontaneously generate one of the other modes (ref. 5). 
(The weak transference of turbulence into noise described 
by the Ldghthill theory is a higher-order effect (ref. 1).) 
The presence of a shock wave, however, provides a mecha- 
nism for a very strong transference: thus, when any one of 
the three modes — turbulence, entropy spottiness, or noise — 
encounters a shock, the interaction will give rise to all three 
modes, in comparable strength, downstream of the shock 
(refs. 2, 4, and 6). 

The first of these cases, shock-turbulence interaction, has 
been investigated at the NACA Lewis laboratory as an 
outgrowth of reference 2 and is reported herein. The anal- 
ysis of the earlier paper was concerned with a single spec- 
trum wave of a turbulent field and was primarily a study 
in gas dynamics. The present paper reformulates the re- 
sults and incorporates them in a spectral analysis; from the 
analysis come the quantitative effects of the interaction of a 
convected homogeneous field of turbulence with an extended 
plane shock front. (Some results of this work are reported 
in abbreviated form in refs. 7 and 8.) The perturbation 
velocity, pressure, temperature, and density distributions 
behind the shock are described in terms of formulas for 
their spectra, correlations, and mean-square values; these 
are separated into the respective contributions of turbulence, 
entropy spottiness, and noise. 

Numerical calculations are presented for the root-mean- 
square values of the pressure (noise) and components of the 
temperature and velocity perturbations for the Mach num- 
ber range of 1 to <*> • one set of calculations refers to iso- 
tropic initial turbulence, another set to strongly axisymmetric 
initial turbulence (lateral perturbations/longitudinal pertur- 
bations » 36/1). The noise pressure level is also presented 
on an acoustic scale for several levels of initial turbulence. 

SHOCK INTERACTION OF SINGLE SHEAR WAVE 

QUALITATIVE DISCUSSION 

According to the Fourier integral theorem, a turbulent 
velocity field can be represented as a superposition or spec- 
trum of elementary waves. A single spectrum wave can be 
interpreted physically as a plane sinusoidal wave of shear- 
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ing motion (e. g., ref. 9) ; a portion of such a wave is shown 
in perspective in sketch (a): 



(a) Wave of shearing motion. 


A similar wave encountering a shock is shown schematically 
in sketch (b ), 



(b) Convection of shear wave through shock: original unsteady-flow 

problem. 

the wave and the shock being viewed “edge-on.” The 
wave-shock interaction is analyzed in reference 2, and 
what follows first is a brief physical account of the main 
results. The wave is supposed to be convected downstream 
by the mainstream with velocity so that it passes through 
the shock. The passage is evidently an unsteady process, 
since the intercepts of the inclined lines — the planes of 
constant phase or wave fronts — move downward along the 
shock; it can be shown that a sinusoidal disturbance ripple 
will move along the shock with the same speed V. 

The unsteady-flow problem may be treated directly (ref. 
4), or it may be converted to an equivalent steady-flow 
problem by superposing an upward velocity V (ref. 2). 
The conversion is illustrated in sketch (c): 


Shear 

and 

entropy 



(c) Transformation to steady-flow problem by superposition of 
velocity V. 


The cross velocity V therein has been chosen so that the 
resultant stream velocity is parallel to the wave fronts in 
the shear wave; the observer then sees what appears to bo 
a steady sinusoidal shear flow passing through an oblique 
shock. This may be called the equivalent oblique shook. 
(Addition of the upward velocity V is, of course, equivalent 
to transforming to a moving frame of reference.) 

Downstream of the shock, the resultant stream flow is 
deflected according to the laws for oblique shocks; the stream- 
lines are the upper lines in the sketch. The vorticity of the 
initial shear wave is convected along these streamlines 
together with the additional vorticity generated by the 
shock. The net result is a refracted, amplified shear wave 
alined with these streamlines. The angle of refraction is 
just the angle of flow deflection of the oblique shock. 

Superposed on the refracted shear wave is an entropy 
wave of the same inclination and wave length. This wave 
arises from the convection of entropy perturbations gener- 
ated at the shock, precisely as the shear wave results from 
the convection of vorticity. The entropy wave is manifested 
physically as a spatial variation of temperature and density 
at constant pressure, by virtue of the equation of state. 

The nonuniform velocity in the shear flow results in a 
nonuniform pressure jump across the shock. The ultimate 
effect is that the shock front develops ripples, modifying 
the pressure variations, and the resultant pressure variations 
propagate downstream as a plane sinusoidal wave (lower 
lines in sketch (c)). 

The character of this wave depends on whether the 
resultant velocity W behind the equivalent oblique shock is 
subsonic or supersonic; this in turn depends on the initial 
wave inclination through V. When W is supersonic, the 
pressure wave is a plane sinusoidal sound wave; it appears 
as a stationary Mach wave pattern in the steady-flow refer- 
ence frame. When W is subsonic, it may be shown that the 
pressure wave, while still plane, is not a simple sound wave, 
but rather attenuates exponentially with distance down- 
stream of the shock; the resultant disturbance velocity is 
not normal to the wave front, and the wave propagates 
relative to the surrounding fluid at less than sonic speed. 

QUANTITATIVE DISCUSSION 

Elementary wave. — Thus far the waves have been dis- 
cussed only qualitatively. Elementary spectrum waves of 
this sort may be expressed quantitatively in the form 

da=dZ a e m (1) 

(All symbols are defined in appendix A.) The wave-number 
vector k is directed normal to the wave fronts and its magni- 
tude equals 2-7r/wave length. The wave amplitude is given 
by the complex quantity dZ a . When a stands for tempera- 
ture, pressure, density, or entropy, these are simple scalar 
waves. When a stands for the components u, v, w of the 
velocity, these are vector waves; two cases may then bo 
distinguished: the waves are either irrotational and com- 
pressible (sound waves) or rotational and incompressible 
(vorticity waves). (See, e. g., ref. 10.) In the first case 
the irrotationahty condition curl a=0 requires that tho 
velocity s and wave vector k be parallel ( u , v , w proportional 
to ki, k 2 , fc 3 , respectively); the sound waves are thus longi- 
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tudinal. In the second case the incompressibility condition 
diy a=0 requires that the velocity a and the wave vector k 
bo perpendicular; that is, 

k l u J rk 2 v J rk 3 w=0 (2) 

Thus, the vorticity waves are transverse and have the 
character of a shearing motion (see sketch (a)); in the dis- 
cussion they have been referred to as “shear waves.” 

Geometric reexamination of prior results. — The shock- 
interaction process for a single shear wave is given quanti- 
tatively in reference 2, but the results are formulated in two 
dimensions. It will be necessary to reexamine the problem 
geometrically in order that the results may be reexpressed 
in three dimensions. 

A perspective view of the initial shear wave in the new 
Xi, Xj, ^ 3 -coordinate system is shown in figure 1. The portion 
of the shear wave shown is on the downstream side of the 
shock front, which is identified with the x- 2) a- 3 -plane. A plane 
passed through the a^-axis perpendicular to the wave fronts 
cuts the shock in the line Or. At a given instant of time 
this Xi, r-plane corresponds precisely to what is called the 
x, y-plane in reference 2. The angle <p of the x u r-plane with 
the horizontal is then the third coordinate in a system of 
cylindrical coordinates. 

In reference 2 the time was eliminated from the equations 
by employing a frame of reference moving with a velocity V 
downward along the shock front, the so-called steady-flow 
frame of reference. In the present paper all results refer 
to a definite instant of time, t=0. Thus, motion of the 
reference frame plays no part, and the results of the earlier 
paper carry over to the present coordinate system on simply 



Figure 1 . — Perspective view of shear wave in relation to reference 

frame. 



Figure 2. — Projective view of shear wave in relation to reference frame. 


replacing x,y by Xi ,r, respectively. The results of the 
transformation are given in the following sections with the 
disturbances reexpressed in nondimensional form according 
to the scheme 

u, v, w=components of velocity perturbation/critical speed 
of sound a* 

p=pressure perturbation/mean static pressure 
p= density perturbation /mean density 
T=temparature perturbation/mean temperature 
In addition, there are other minor respects in which the 
notation has been modified from that of reference 2; for 
example, the waves are expressed in complex form. 

Initial shear wave (~ initial turbulence). — At time t= 0, 
the velocity field of the initial shear wave is, in cylindrical 
coordinates, 


du=dZ u e ' 
dv T =dZ T e i * i ► 
dv r =dZ f e i ^- „ 


( 3 ) 


where du is parallel to x 1 (longitudinal direction), dv T is parallel 
to r, and dv 9 is perpendicular to r and Xi, in the direction of 
increasing <p (see figs. 1 and 2). The wave-number vec- 
tor k lies in the x h r-plane, making an angle 0 with the r- 
axis. 

Refracted shear-entropy wave (~ final turbulence and 
entropy spottdness). — The velocity field of the refracted 
shear wave (fig. 3) is 


du'=dZ' u e 

dZ' u —X dZ n 

dv' r =dZ‘, 

dZ'=Y dZ r 

do' r =dZ'^ 

dZ' r =dZ 9 


( 4 ) 


at time f= 0, where k' is the new wave-number vector, making 
an angle 0' with the r-axis. The radial components of k' 
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and k are equal (k' r =k T ), and the further dependence of 
k' on k is expressed through the dependence of 0' on 0. 
Similarly, the complex amplification factors X and 7 de- 
pend on k in terms of 0. Expressions for X, 7, and 0' are 
given in appendix A* 

The perturbation pressure dp' will be zero because this is 
again a shear wave, free of accelerations. The temperature 
perturbation associated with the companion entropy wave 
(fig. 3) will be 


dr'=dZ f T e i ^ r '- dZ' r =T dZ u (5) 


With p '— 0 (to the first order), the dimensionless density 
perturbation p' will be just the negative of the dimensionless 
temperature perturbation r', according to the linearized 
equation of state. The form of the function X is given in 
appendix A. 

Aside from the change in wave inclination, the description 
of the refracted shear-entropy wave in terms of the initial 
shear wave depends entirely on the amplification factors X 
and 7 and the function T. Such functions play a role 
similar to the “transfer functions” of the theory of servo- 
mechanisms (ref. 11), and it appears appropriate to carry 
the name over to the present field. 

Generated sound wave (~ noise field). — The shear-entropy 
wave downstream of the shock is accompanied by a plane 
irrotationai pressure wave (sound wave) of different in- 
clination (see fig. 3). For small inclinations 0 of the initial 
shear wave, this pressure wave attenuates exponentially with 
distance from the shock; for inclinations greater than a 
certain critical value 6 cr (see appendix A), the pressure wave 
is unattenuated. The critical wave inclination 0 cr corre- 
sponds to the attainment of sonic speed in the mean flow 
behind the “equivalent oblique shock” referred to in the 
qualitative discussion. 

The velocity field can be represented in the form 


du"=dZLe^ 

dvl^dZU*?-*- 

<K=dzr 9 e^ 


d2Z=xdZ u " > 
dZT r =TdZ u j- 

dz;=o 


( 6 ) 


where k " is the wave-number vector, making an angle 0" 
with the r-axis; again the radial component matches that of 
k ; namely, k"=k T . The sound-wave angle 0" and the 
transfer functions x and X are specified functions of the 
shear-wave angle 0; moreover, for O<0<0cr, X and X are 
functions of x u showing an exponential decay to zero as 

Zr- 

The pressure perturbation may be written 


dp"=dZle^ dZ' v =PdZ u (7) 


where P=P(zx) is a transfer function defined in appendix A; 
like x and X, P decays exponentially with x for O<0<0 cr . 
The corresponding density and temperature perturbations 
are proportional to p"; they may be obtained from the 
isentropic property of the sound wave as p"=p'77 and 
T r, =p"(y l)/x- 


Transformation to Cartesian coordinates. — Expressions 
for the velocity field in Cartesian coordinates will be needed. 
The transformation from cylindrical coordinates is effected 
by means of the relations 

dZ B =dZ r cos <p—dZ 9 sin <p 
dZw—dZ T sin (p-^-dZ# cos <p 



where primes (~ refracted shear wave) or double primes 
(~ sound wave) may be inserted throughout as needed. 

The transformation results in 


Initial f 
shear ^ 


dct=dZ a e i ±'- ) where a=u, v, w 


} 


(9) 


C da' =dZ' a e i *''-, where the values of dZa^\ 
for a=u, v , w are, respectively, 


Final dZ u — XdZ u 
Shear 

dZ' 9 =7 dZ r cos <p—dZj sin <p 


y do) 


L. dZ' w —7 dZ T sin (p+dZ? cos <p 


f da"=dZa6 i ?-z, where the values of dZ* a ' 
for a=u, v, w are, respectively, 


Sounds 


dZ:=xdZ u 
dZ %= T dZ T cos <p 


Y (ii) 


[^dZl^TdZr sin (p 

SPECTRAL ANALYSIS OP RANDOM FIELDS 


The foregoing relations will be fitted later into a spectral 
analysis of the fields of turbulence and noise. Appropriate 
analytical techniques can be found in the spectral theory 
of random functions; suitable developments of this sort are 
given by, for example, Moyal (ref. 10) and Batchelor (ref. 
12) for spatially homogeneous fields. The first part of the 
present section will be devoted to an interpretation (with 
some liberties) of relevant parts of the two papers; the latter 
part will be devoted to developments for inhomogeneous 
fields and for correlations of a two-dimensional field with 
a three-dimensional field. 



Figube 3. — Interaction of shear wave with shook: view in X\ t r-plano, 




SH0CX-TTJRBT7LBNCE INTERACTION 
HOMOGENEOUS FIELDS 

Amplitude spectra. — Consider a three-dimensional field 
of small disturbance (e. g., turbulence or noise) of unlimited 
extent. Let this field be homogeneous in the sense that the 
statistical properties do not vary from point to point. The 
instantaneous spatial distribution of any physical quantity a 
can then be represented mathematically by a Fourier integral 
in the Stieltjes form (refs. 10 and 12) 

a(x)=j e^'-dZalk) (12) 

where the triple integral goes from — ® to ® in each com- 
ponent of k= (k h k 3 , k 3 ). 

If equation (12) is written in the form 



then, by comparison with equation (1), da can be identified 
with what has been called an elementary spectrum wave. 
The Fourier integral is thus to be interpreted as a superpo- 
sition of infinitely many of such plane waves. In the inte- 
gral the components of k take on all values; it follows from 
the significance of k as a wave-number vector that all wave 
inclinations and wave lengths appear. An aggregate of vor- 
ticity waves with a suitable distribution of amplitudes among 
the various wave lengths and inclinations can represent a 
turbulent field (ref. 13). Similarly, an aggregate of sound 
waves suitably distributed can represent a random noise 
field (ref. 10). Finally, an aggregate of the scalar entropy 
waves can represent a random field of entropy spottiness. 
A combination of these three basic types of disturbance — 
entropy spottiness, turbulence, and noise — constitutes the 
most general random small-disturbance field that may exist 
in a gas (refs. 4 and 5). 

Correlations. — Let a be measured at some point P and /3 
at some point P a vector distance £ from P; then the space 
average of the product a/3 as P and P vary but their vector 
separation £ is held fixed may be defined as the space-average 
correlation a/9(£). Alternatively, the disturbance field may 
be considered to be just one of a large number, or ensemble, 
of statistically similar fields (e. g., the flow fields of a great 
many “identical” wind tunnels operated simultaneously); 
the average of a/3, with P and P fixed, over all members of 
the group, is the ensemble-average correlation. The equa- 
tions that follow, from the theory of random functions, refer 
solely to ensemble averages, but space averages are desired 
in practical applications. The ergodic hypothesis of proba- 
bility theory equates the space average to the ensemble 
average provided that, at any instant, the disturbance fields 
a and /3 are stationary random functions of position; that is, 
the disturbance fields are spatially homogeneous. 

In what follows, the term “cross-correlation” will be ap- 
plied for a 3, the term “self-correlation,” or simply 
“correlation,” for a=/3. 

Correlation and power spectra. — The cross-correlation 
a/3 (£) (like a or /S, indvidually: see eq. (12)) may be ex- 
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pressed by means of the Fourier integral as a spectrum of 
plane sine waves: 

^3(|)=Je«-l[a/3 ]dk (13) 

where [ap] is a function of k, and dk is an abbreviation for 
dkidk^dkz. The differential e^*[ap]dk may be regarded as 
the contribution to the correlation made by spectrum com- 
ponents with wave number between k and k+dk . The 
function [aP] is called the “spectral density” when a=p, the 
“cross-spectral density” when a^P (ref. 11). The array of 
nine spectral densities signified by [aP] when a and P are 
limited to mean u, v, or w is the “spectral tensor” of the 
velocity field and is commonly written as of (The 
corresponding array of the nine velocity correlations aP (£) 
is the “correlation tensor,” commonly written as i? y (£).) 

Equation (13) includes as a special case the self-correlation 
or mean-square relation 

\aa]dk , where £ — 0 (14) 

If a were a velocity component (say u), then a? would be 
twice the space-average kinetic energy per unit mass associ- 
ated with u. The spectral density [act] is in this case an 
energy density (per unit mass, per unit wave number). For 
similar reasons, where spectra of the kind defined by equation 
(14) have occurred in physics (e. g., in the harmonic analysis 
of radio noise), they have generally been called energy, 
intensity, or power spectra. 

Correlation spectrum in terms of amplitude spectra. — 
The rather analogous forms of equations (12) and (13) are 
of interest. Equation (12) expresses the spectrum of the 
amplitude of the fluctuating quantity a; this may be termed 
an amplitude spectrum. Equation (13) expresses the 
spectrum of the correlation of a with p; this has been termed 
a correlation spectrum. The complex magnitude dZ a Qc) of 
the amplitude spectrum fluctuates in an apparently random 
manner as k is varied (refs. 10 and 12). The magnitude 
[ap\dk of the correlation spectrum, on the other hand, varies 
smoothly with &, since the correlation is a smoothed or 
averaged quantity (ref. 12). The amplitude spectrum gives 
no direct information concerning averaged (i. e., statistical) 
properties of the disturbance field, whereas the correlation 
spectrum leads directly to expressions for correlations and 
mean-square values (see eqs. (3) and (4)). One-dimensional 
spectra and scales of turbulence can also be determined 
(e. g., ref. 14). 

It would be desirable to formulate the shock-turbulence 
interaction problem directly in terms of correlation spectra, 
but formidable difficulties stand in the way. It has been 
simpler to start with the shock interaction of a single shear 
wave, which deals with amplitude spectra, and to infer from 
this the changes in the correlation spectra. The whole 
procedure depends on the following relation (refs. 10 and 12) 
which connects the two kinds of spectra, namely, 


[ap]dk=dZ*(k)dZp{k) 


(15) 
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where dZ(Jc) is associated with the wave-number range 
between k and and the bar represents the ensemble 

average. This relation is fundamental to the spectrum 
analysis of the present paper. Its significance is this: the 
single-wave analysis (summarized in an earlier section) pro- 
vided the change in amplitude of an individual spectrum 
wave in the form dZ a -*dZ' a , say, and similarly, dZp-*dZ' f ) 
equation (15) provides the means for determining therefrom 
the corresponding change in the spectral density: [a0]-»[a'/?']. 


random functions. The simplification can still be achieved 
by replacing the inhomogeneous fields by “equivalent” 
homogeneous fields that match, respectively, at the points 
x and x'. This is accomplished by freezing dZt(k,x) in 
equation (17) at the value dZ*(k,x) while allowing x to vary 
in the exponential, and correspondingly freezing dZp in 
equation (18). 

When applied to the so-defined equivalent homogeneous 
fields, equation (20) reads 


INHOMOGENEOUS FIELDS 

The spectral representation of a spatially homogeneous 
random field is given by equation (12) : 

a(x)=j e**dZ a (k) 

A corresponding possible representation of an inhomogeneous 
field is 

a(x)=J e^dZaQcjX) (16) 

where dZ a now depends on position; the gound field behind 
the shock is of this character. The following spectral 
analysis of such inhomogeneous fields is a development of 
Moyal’s treatment of homogeneous fields (ref. 10). 

Let a(x) and fi(x') be inhomogeneous fields 


a(x)= 

(17) 

PW-Je^dZ^xT) 

(18) 


where equation (17) is an alternate form of equation (16). 
The correlation of a and /3 for fixed positions x and x', re- 
spectively, can be formed by taking the ensemble average 
of their product: 


«(£) /3&0= JJ e* dZ*(k,x) dZ^Qc', x') (19) 


dZXkjX) d^(X- / ,iO= [$1 dkcUd 6 (k'-k) (21) 

where the ^ over [ap] signifies the functional dependence on 
x and x'. Upon substitution into equation (19) and integra- 
tion over V there results, with \ 

eti&pfg) =J e*-l M dk (22) 

The spectral density [«/3] can be evaluated by integrating 
equation (21) over V : 

[$] dk=dZ* (k, ®) dZp (k, i') (23) 

where the integral property of the 5-function, 

J_" f(k r )S(k!-Jc)dk!=M 

has been used, with/(&') an arbitrary function. 

Equations (22) and (23) for inhomogeneous fields are of the 
same form as their counterparts, equations (13) and (15), 
respectively, for homogeneous fields. In the homogeneous 
case the dZ’ s are functions of position, and equation (23) 
mplies a corresponding dependenc e of [ap] on position. 
Moreover, the correlation a(x) P(& 0 depends on | and & 
iseparately as well as on their separation 

CORRELATION OF TWO-DIMENSIONAL FIELD WITH THREE- 
DIMENSIONAL FIELD 


The operations of integration and averaging commute, so the 
averaging bar may be regarded as placed over the dZ ? s alone 
on the right side. 

Equation (19) could immediately be simplified if the fields 
a(x) and p{xf) were homogeneous; in that case the important 
relation 


The local perturbations of the shock face from the moan 
(x£, x£) plane constitute a homogeneous two-dimensiona 
field of the general form 

P(x’ u 4, 4)=f e i( ^ +iitS> dW fl (4, *S) (24) 


where 


and 


dZt(@dZ ff Qd)=[afi dkdk' 5(k'-k) (20) 


8(k?—k )= 0 iovld^k 
= co for k'=k 

f 0 8Qd-k)dk '= 1 

— co 


would hold (ref. 10), according to the spectral theory of 


where x[ has been fixed at the value It may be desired 
to correlate such a field locally with a three-dimensional field 
(e. g., the turbulent velocity field). To this end, equation 
(24) is rewritten in the form 

dW4KX) 

Now, if xi in e ~ <fc # is fixed at the value x[ } p will be general- 
ized to a three-dimensional field (elementary wave number 
Jc') that matches the original two-dimensional field in its 
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plane of definition x 1 =x 1 . This “equivalent” field may be 
written 


where 


/ 9(*0 =Je a -’- x -' dZtQcT) 

dZg(kf)=e~ tt '^ dWM, ^ , 


(25) 


Equation (25) is of the form of a three-dimensional homo- 
geneous field and may be used in place of (24) in equations 
(13) and (15) to provide the correlation of 0 with any three- 
dimensional homogeneous field in the common plane Zi—Zi. 


and also, by the coordinate transformation (8), 
dZ*dZ r +dZ*dZ^=dZ*dZ B +dZtdZ w 
Thus, equation (27) becomes 

dZ?dZ:+dZ'SdZ' a =(\r\ 3 -l) tan 3 Q dZ&Z u +dZ*dZ 9 +d2S#Z v 

(28) 

Application of equation (15) yields 

[u'u'] dtf=\X\*[uu]dk \ 

{ b V] + [w f w r \ } dfc— (\Y\*~ 1) tan 3 6 [uu] dk+ { [vv] + [ww] } dk j 

(29) 


SHOCK INTERACTION OF SPECTRUM OF SHEAR WAVES 
(TURBULENCE) 

The interaction of a single shear wave with a shock has 
been discussed in detail. With this as the basis, the statistical 
behavior of a spectrum of shear waves representing turbu- 
lence will now be derived; the procedure will make use of 
the spectral-analysis relations of the last section. The 
problem is formulated as follows: given the spectra (and 
hence correlations and mean-square values) associated with 
the turbulence convected into the shock, to calculate there- 
from the spectra, correlations, and mean-square values 
associated with the turbulence, entropy spottiness, and 
noise in the flow downstream of the shock. 


These are the desired expressions relating diagonal elements 
of the spectrum tensors of the turbulence on opposite sides 
of the shock. 

Noise field. — If operations similar to those of the last 
section are applied to equations (11), there results 


[u''u'^dV f =\x\ % [uu]dk 1 

{[t/V'l +[^'W / ]}d^ , =|T| 3 tan 2 ^M(a J 


(30) 


These equations relate the diagonal elements of the spectrum 
tensor of the noise generated behind the shock to the longi- 
tudinal spectral density of the initial turbulence ahead of 
the shock. 3 


DIAGONAL TERMS OF VELOCITY SPECTRUM TENSOR 


MEAN-SQUARE VELOCITY COMPONENTS 


The respective spectrum tensors for the turbulence and 
noise downstream of the shock each consist of nine elements; 
of these the three diagonal terms are most important since 
they lead to the mean squares of the velocity components. 
The relatively simple derivation of the first diagonal term 
and the sum of the second and third will be carried out in 
the present section. The derivation of the complete tensor 
is carried out in appendix B by a more formal procedure. 

Turbulence field. — The shock interaction effects have 
been expressed in terms of relations between wave amplitudes 
on opposite sides of the shock (eqs. (9) and (10)). Cor- 
responding relations between spectral densities (elements) 
on the two sides can be obtained by use of equation (15). 
Some preliminary manipulation is required; thus multiply 
both sides of equations (10) by their complex conjugates, 
and add the last two; there results 

dZ' u *dZ' u =\X\*dZ*dZ u (26) 

dZ'*dZ'+dZ' w *dZ' w = \Y\ 2 dZ r dZ*+dZ*dZ p (27) 

But by geometry (fig. 2), 

dZ T dZ u tan 6 


Turbulence field. — The mean-square velocity components 
follow directly from integration of the spectral density (see 
eq. (14)). Integration of both sides of equations (29) yields 


® ,s 4-'to ,s =?+'W? s +J' (|F| J — 1) tan 2 0 [««]<& 


(31) 


Thus, the mean-square velocity components behind the 
shock (primed values) are given in terms of those ahead of 
the shock, the single-wave transfer functions X and Y, and 
the longitudinal spectral density [uu] of the initial turbulence. 
Note that X and Y are functions of k in terms of 6 (see 
appendix A). 

Noise field. — Similarly, integration of equations (30) 
yields the mean-square velocity components in the noise 
field: 


[x | 2 [uu]dk 

v FT *-\ r v/' 2 = J |T| 2 tan 2 d[uu] dk 


(32) 


dZf=dZ* tan 0 


Here again, x and T are functions of k in terms of 0. 


t Direct expressions for the spectra downstream of the shock may be desired, free of the unequal volume elements dk f dkf, or dk". This may he effected in eq. (29) by dividing both sides 
by dkf; then (since dk Is shorthand for dkidktdk i, and similarly for dk') the ratio dkjdV may be interpreted as the Jacobian (say JO for the transformation from ktotf. Upon evaluation, 



Similarly, In eq. (30) divide by dk" and interpret dkfdk " as the Jacobian (say J") for the transformation from k to Upon evaluation 


J"- 


cos* $" 
m cos 1 O' da" 
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MEAN-SQUARE PRESSURE 

The first-order pressure field is associated solely with the 
noise field: the pressure field associated with the turbulence 
is of the second order in velocity and may be neglected in 
comparison. 3 The spectral density of the noise pressure 
can be related to the spectral density of the longitudinal 
velocity in the initial turbulence; the relation is obtained by 
multiplying both sides of the second of equations (7) by 
their complex conjugates, averaging, and applying equation 
(15) to each side: 

kp"p"] \P\ 2 [uu]dJc (33) 

The integration of both sides of equation (33) yields the 
mean-square pressure in the noise field as 

^h=J| p\*[uu]dk (34) 

MEAN-SQUARE TEMPERATURE 

The temperature perturbations in the noise field, because 
of the isen tropic relation, are equal to (y— 1 )/y times the 
pressure perturbations; thus, the relations corresponding to 
equations (33) and (34) may be written down at once. 

The temperature perturbations associated with the entropy 
spottiness behind the shock require a separate analysis. 
The spectral density of the temperature perturbations can 
be evaluated by operating on equation (5) in the now- 
familiar manner (see remarks preceding eq. (33)) ; the result 
is 

[Vt 7 ] dJc , = \T\ 2 [uu] dk (35) 

The integral relation obtained from equation (35) is 

7*=f (36) 

This equation evaluates, for the region behind the shock, 
that part of the mean-square temperature spottiness asso- 
ciated with the entropy spottiness. 

MEAN-SQUARE DENSITY 

It is unnecessary to write down special expressions for the 
density field: the respective contributions of entropy spotti- 
ness and noise to the density perturbations are related to 
the corresponding temperature and pressure perturbations 
by p'=— r' and p"/ 7 , according to the small-perturbation 
form of the equation of state. 

CORRELATIONS NOT JOINTLY INVOLVING TURBULENCE AND NOISE 

Attempts at simplification. — If the spectral density 
[o0](£) is known, the corresponding two-point correlation 
a/3 (£) can, in principle, be obtained by means of equation 


(13). In this fashion, for example, the longitudinal velocity 
correlation in the turbulence behind the shock may be 


expressed, with use of equation (29), as 


ttV(£0 =J \X\*[uu]eV-m 

(37) 

or 


=J|X|*[tm] «**'•!' J'dk' 

(38) 


(See footnote 2, p. 7, for significance of J f .) 

Either of the forms (37) or (38) may prove awkward 
because of the admixture of Jc and k' in the integrand (e. g., 
[uu] is ordinarily most simply expressed as a function of Jc ) . 
However, it is possible to find a fixed vector £ that satisfies 
the relation &' •£'=£*£; this gives the more convenient 
relation 

uV(£'j=f |X| S [uu] e*idk (3 9) 

where £=m£, &={£, &=£. In all the self- or cross- 
correlations involving properties of the turbulence and 
entropy spottiness behind the shock, whether they be 
velocity components, temperature, density, or entropy, the 
transformation &' .£'=£*£ can be made to simplify the 
exponential. 

The physical- interpretation of the relation between j? and 
£ 7 is this: if two fluid particles upstream of the shock are a 
vector distance £ apart, after convection through the shock 
they will be a vector distance £ 7 apart. Put another way, a 
“box” of turbulent fluid of edges &, £3 will be compressed 

on passing through the shock and will emerge downstream as 
a shorter box of edges £. Therefore, equation (39) in 
effect expresses correlations in the space downstream of the 
shock in terms of equivalent correlations in a stretched space 
upstream of the shock. 

The analog of equation (37) for the correlations of proper- 
ties of the noise field involves £ // *£ // in the exponential, 
rather than Je'-jt'. Here no great simplification appears to 
be possible in general: 4 there exists no fixed vector £ that 
satisfies the relation &".£ ,/ =&.|. This lack reflects the 
nature of the transformation from Jc to Jc": the respective 
components of the two vectors are not in fixed proportions, 
but instead vary with the inclination of Jc. The particular 
coordinate compression £— that works for the turbulent 
field (it expresses the change in dimensions of a fluid “box” 
convected through the shock) will not work for the noise 
field. An exception occurs when £ 7/ is chosen parallel to the 
shock plane (radial direction, ^=0). Then 
and since Jc^ ~k rJ it follows that for this case & 7/ .£" = hg* . 

The integral for a particular correlation simplifies con- 
siderably when £ (or £ 7 , or £") is taken in the direction of one 


J The local pressure field associated with turbulence, although weak by aerodynamic standards, may be strong by acoustic standards. If the turbulence (e. g., In a boundary layer) Is 
convected past a stationary microphone, a strong response can be observed; the phenomenon Is called "psuedo-eoundL" The noise sensation produced by wind blowing past the cars Is 
presumably a similar effect associated with turbulent separation of the flow. 

’ * A partial simplification Is k-V'+iklf-kdiY. 
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of tlie coordinate axes, say z t . In the former case, b £ 
becomes £<£; and the exponential can be replaced by cos 
k&t since the imaginary sine component will integrate out. 
Similarly e can be replaced by cos k 

Cross-correlations. — The phase angles of the transfer 
functions must be considered in formulating cross-correlations. 
For example, the correlation of local temperature with 
longitudinal velocity in the entropy and turbulence fields 
behind the shock is readily obtained as 

Fw(0 )= J {Te {l r)*{Xe il ')luu]dk 
= J TXe i(l >-^ [uu] dk 

The integrand, except for the exponential, is even in the 
wave inclination 6; the phase angles 8 , and 8 T (in the notation 
used) are odd in d (both properties can be inferred from the 
symmetry of the wave-refraction process with respect to 6 ). 
Accordingly, the imaginary sine term in the exponential will 
integrate out, and 

r 7 iZ 7 (0 )= J TX cos (8 B —8 T )[uu]dk (40) 

The corresponding relations for other cross-correlations can 
be written down by analogy. 

CORRELATIONS BETWEEN TURBULENCE AND NOISE 

Cross-correlations between the turbulence and noise fields 
require a special treatment, partly because of the inhomo- 
geneity of the noise field, and partly because of the non- 
parallelism of the physically associated waves. In what 
follows, an expression for the correlation of noise pressure 
with longitudinal turbulent velocity will be derived. From 
this the qualitative variation of the correlation with distance 
downstream of the shock will be inferred. 

The refracted shear wave (~jfc') and pressure wave (~k") 
associated in an elementary interaction process have different 
inclinations (fig. 3). As a consequence, the formal applica- 
tion of the relations given in the section SPECTRAL 
ANALYSIS OF RANDOM FIELDS leads to difficulty: the 
spectral density of any correlation appears to vanish accord- 
ing to equation (21). Actually, the formulas are inapplica- 
ble to correlations involving mutually inclined waves; this 
will be brought out clearly in the following derivation of the 
applicable formulas. For simplicity the derivation will be 
limited to the correlation of turbulent longitudinal velocity 
v/ at point x' with noise pressure p" at point x"; extensions 
to other cases are straightforward. The derivation will first 
be carried out as though the noise field were homogeneous 
(no variation of transfer function P with x) } and then will 
be adapted to take account of the actual inhomogeneity. 

The respective Fourier integrals may be written 

*'(£)- J e t ^'dZ u (k')=J, e-v-z’dZtQd) 


P"<2?)= j e‘ k ' x -dZ v {P') 

The correlation may be formed as the ensemble average of 
the product u'p" : 

u , (^)p"(^)=jj ~^dZ*(k0dZ p (t) (41) 

where the bar has been taken inside the integral, since the 
operations of averaging and integration commute. Equation 
(7) and the first of equations (4) may be used to simplify 
the right side: 

dZ*Qd)dZ v (fc) =X*Qc)P(P)dZi(}c)dZ u (P) (42) 

A A 

■where k " bears the same relation to & as k " does to k. By 
virtue of equation (20), equation (42) reduces further to 


dZ*(Jd)dZ r (t)=X*(]c)P(i) [uu]5(i~k)dk 4 
if the fields are homogeneous. Substitution of this relation 

A 

into equation (41) and integration over k result in 

« / (x0p"(2 / 0=J e* (43) 

A 

since the 5-function eliminates all values of k but k and 

A 

similarly all values of k" but k". Finally, the equation may 
be generalized to apply to the actual inhomogeneous pressure 
field, according to equation (23) and the discussion preceding 
it, by writing P (Jfc) as P (k, x') and using the value appro- 
priate to x[. 

Equation (43) is the general relation for the two-point 
correlation of longitudinal turbulent velocity u r with noise 
pressure p n . The striking feature is the difference of the 
exponential term from those in equations (13) and (22) ; this 
constitutes an a posteriori demonstration of the inapplica- 
bility of those equations. 6 j 

If the turbulent velocity and noise pressure are correlated 
locally (z"=x'), the expression simplifies to 

W T (^)=J e^'^XQcjPQc, aft [uu] dk (44) 

since Directly at the shock, and the 

right side simplifies further; the integration can readily be 
carried out for isotropic turbulence, and a nonvanishing 
correlation will be obtained. Behind the shock (a^>0), 
the exponential oscillates sinusoidally; for a given wave 
inclination the behavior is essentially like cos Ckx[ y where 
G is a constant. For x[ very small, the cosine is near unity 
over the significant range of k (the range for which [uw]^>0). 
Hence the correlation is only slightly diminished at small 
distances behind the shock. At somewhat greater distances 
the oscillatory nature of the cosine begins to be felt before 
[uu] dies out, and the correlation falls off noticeably. Finally, 
at very large distances, cos Okx[ oscillates over a great 


* However, eq. (43) la equivalent to that which would result from eq. (13) or (22) upon replacing the pressure wave by a locally equivalent shear wave parallel to the actual shear wave, os 
discussed In ref. 8. 



692 


REPORT 1233 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


many periods as k covers its important range, and the plus 
and minus contributions to the integral cancel each other; 
thus at these large distances behind the shock the noise- 
turbulence correlation falls to zero. 

INTERACTION OF TURBULENCE WITH AN OBLIQUE SHOCK 

All the foregoing analysis may be applied to an oblique 
shock by treating the latter as a normal shock with a super- 
posed cross-velocity which is to be ignored. The coordinate 
system should be oriented so that the 3 i-axis is normal to 
the oblique-shock front (on the downstream side), and the 
av- and Tj-axes lie in the ^ shock front with the 23 -axis in the 
plane of the stream-velocity vector and the sraxis. The 
component of the stream velocity in the 21 -direction is the 
U velocity of the equivalent normal shock. From here on 
the analysis for the normal-shock case may be applied. 

Ordinarily the turbulence spectrum tensor will be defined 
(as & i0 say) in a system x[, a 4 , with the a^-axis alined 
with stream direction, and it will be necessary to transform 
& if to the new system x u 23 , 23 . If the shock angle of the 
oblique shock is the primed and unprimed axes are related 
according to the following scheme: 



Xi 


XJ 

*5 

r u— sin £ 

ru— 0 

m—cos yf' 



Tn=l 

TO-0 

j' 

i 

ni“—cos f 

m— 0 

m— sin $ 


The corresponding element of equation (48) is 

w iw+m (si) 

Thus the tensor elements $ n and have the same functional 
form, reflecting the isotropic property of invariance under 
rotation of coordinates. This particular example of the 
coordinate rotation applied to isotropic turbulence is trivial 
in that the result could have been written down in advance 
without recourse to the transformation equation. Never- 
theless, it illustrates the formal application of the trans- 
formation and, in addition, serves as a check on the first of 
equations (47) in yielding the required invariance, 

CALCULATIONS 

Numerical calculations have been carried out for flows 
in which the turbulence incident on the shock iB (1) isotropic 
and (2) strongly axisymmetric. An account of the iso- 
tropic case follows. The more complicated axisymmetric 
case adds little of interest and is therefore left to appendix 0. 

MEAN-SQUARE VELOCITY COMPONENTS IN TURBULENCE FIELD 

The equations that jointly relate the upstream (unprimed) 
and downstream (primed) mean squares ore 

[uu] dk (52) 

(53) 


where is the cosine of the angle between x[ and x The 
transformation is effected by the formula 


I'M-to'*- f |S| 2 sin * 6 - [uu] dk+S+tf (54) 

J COS o 


4 > M =r in r^4>,'j (46) 


where the repeated indices i and j are to be summed over. 
The diagonal terms in the result axe relatively simple: 


$ u =$a sinV+^M cosV-sin ^ cos tft ($xs+$ji) ") 




Y (47) 


The first of these is just equation (14) with a=u; the last two 
result from substituting into equations (31) the expressions 
for |X]* and |F] S from appendix A. So far the equations 
have not been specialized to isotropic initial turbulence. 

When the initial turbulence is isotropic (i. e., has spherical 
symmetry), its longitudinal spectral density [uu] has the 
general form (e. g., ref. 12, eq. (3.4.12)) 


$ S 3 =$[ 1 cosV+$« sinV+sin $ cos $ ($[j+$[i) J 


[uu]=k?F (k) cos 2 0 


( 66 ) 


The coordinate transformation whereby goes over 
into $ nn may he illustrated most simply by choosing & it to 
correspond to isotropic turbulence; in that case, has the 
general form (e. g., ref. 12) 

& <i =FQc r )(Jc n 6 { ,-klk' i ) (48) 

Substitution into the first of equations (47) yields 

=F(k r ) [(£**+&*) sin 2 *+(&*+&*) cos 2 ^+2 k'M sin * cos <p] 

=F(k') [&*+(£ sin t+Jc[ cos m (49) 

In the preceding equations, k[, k' U! k' t are the components 
of the wave-number vector in the primed coordinate system; 
these are related to the components k i} k a , k z in the unprimed 
system precisely as x[, *£, x z axe related to ®i, x 2 , x z in equa- 
tions (45). As a consequence, equation (49) can be readily 
shown to reduce to 


where F(k) is an arbitrary function of k. ( FQc ) will ulti- 
mately cancel out in forming ratios.) It is appropriate, 
then, to go over to a form of spherical polar coordinates: 

k!=—k sin 6 
k a =k cos 0 cos <p 
k z =k cos 0 sin <p 
dk=k f cos d dk dtp dd) 

Equations (52) and (53) may now he -written 

r™ rsr ri/3 

u 2 =2 k*F (k) dk I (bp I cos 3 0 dd ( 57 ) 

Jo Jo Jo 

pit fin 

'ii /2 =2j k*F(k)dkj dtp\ |<S| 2 cos 2 0' cosOdB (68) 


y ( 56 ) 




(50) 


where the factor of 2 and the limit 7r/2 result from the sym- 
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metry in 6. Division of equation (58) by (57) yields, since 

rr/2 

I cos 3 0 d6= 2/3, 

u /2 /u 5 =|J > |£| 2 cos 2 0'cos 0 d6 (59) 

In a rather similar fashion, equation (36) yields 

tj' 2 * /^==^' 2 /t?==| ^1+J* l^l* 8 ^ 2 ^ 008 Q (60) 

where use has been made of the initial isotropy u 2 =v 1 — / w z i 
and final axisymmetry v ,2 =w' 2 . 

The transfer function S in equations (59) and (60) is a 
measure of the amplification of a single spectral component 
in passing through the shock; the associated phase angle is 
5, (not relevant here). S , like the other transfer functions, is 
a complicated function of 6 that does not lend itself to ana- 
lytic integration. A numerical tabulation of S and 5, against 
0 is given in tables I (c) to (k) for the respective Mach num- 
bers of 1.10, 1.25, 1.5, 2.0, 2.5, 3.0, 4.0, 6.0, and <*>; these 
tables were used in conjunction with numerical integration 
to evaluate equations (57) and (58). (S reduces to 1 for all 
6 at 

MEAN-SQUARE TEMPERATURE IN ENTROPY FIELD 

The derivation of r 75 /?? is parallel to that of v^fu 2 ^ equa- 
tion (53) being replaced by equation (36). The result is 
(analog of eq. (59)): 

?5=| iff''* |T|* cos 8 e de (61) 

The transfer function T and the associated phase angle 
8 r (not relevant here) are tabulated against 0 in tables I (c) 
to (k) for the various Mach numbers. The tabulated values 
used in the numerical integration of equation (61). 

MEAN-SQUARE PRESSURE IN NOISE FIELD 

Because of the similarity of equations (34) and (36), the 
mean-square pressure can be written down by inspection of 
equation (61): 

|P| 3 cos 3 6 dd (62) 

The integration has been performed numerically with use 
of the definition of P in terms of n (appendix A) and the 
values of n against 0 tabulated in tables I (a) to (i), appro- 
priate to x= co . Thus, the integral as evaluated refers to the 
asymptotic mean-square pressure far behind the shock. 

RESULTS AND DISCUSSION 

The results of the calculations of the preceding section are 
show in figure 4 for Mach numbers of 1 to ® ; this figure 
evaluates the disturbance field — both turbulence and noise — 
downstream of a shock when isotropic turbulence is eon- 
vected into the shock. The velocity perturbations, on a 


root-mean-square basis, are in percent of stream velocity 
ahead of the shock (thus the basis is the same on both sides 
of the shock); the temperature and pressure perturbations 
are in percent of ambient behind the shock. 6 The velocity 
curves refer solely to the turbulence component, the tem- 
perature curve to the entropy component, and the pressure 
curve to the noise component of the field behind the shock. 

The curves show that isotropic turbulence is somewhat 
transformed in passing through a shock, the longitudinal 
and lateral components no longer being equal; the selective 
effect is, however, mild compared with that of screens or 
wind-tunnel contractions (cf., e. g., ref. 14). In addition, 
although the incident flow was assumed isothermal and isen- 
tropic, the downstream flow possesses an entropy spottiness, 
which is a “frozen” convected pattern like the turbulence. 
The root-mean-square temperature associated with the en- 
tropy spottiness, in percent of ambient, is seen to be not 
much less than the root-mean-square velocity of the initial 
turbulence, in percent of free stream. 

In the theory the entropy spottiness is spatially correlated 
with the longitudinal component of the turbulent velocity 
everywhere behind the shock. In practice it is to be ex- 
pected that the correlation will soon be destroyed by eddy 
intermixing as the combined fields are convected down- 
stream from the shock; this intermixing, being second order, 
is neglected in the linear theory. Directly at the shock, the 
noise pressure likewise is correlated with the longitudinal 
component of the turbulent velocity. According to the 
earlier qualitative examination, however, this correlation 
falls off with distance behind the shock, reaching zero far 
back. 

The peculiar hump in the curve of root-mean-square 
noise pressure against Mach number just above 1 has 
commanded special attention. In order to delineate the 
shape accurately, additional numerical computations (be- 
yond those for the other curves) were made at Af— 1.05 and 
M =1.01. These were supplemented by an analytical study 
which established that the curve varies like (M— 1) 1/4 as 1 



Figure 4. — Disturbances produced behind shook by interaction with 
isotropic turbulence. Turbulent intensity just before shock, 0.1 
percent. Root-mean-square velooity in percent of initial stream 
velocity ahead of shock; root-mean-square temperature and pressure 
in percent of ambient behind shock. 


4 For the circumstances of figs. A and 5, namely, longitudinal component of initial turbulence equals 0J. percent of stream velocity, the dimensional quantise* plotted are as follows In terms 

of the non dimensional symbols used In the analysis: longitudin al velocit y, percent initial stream velocity, OJ. V tF*/t**; lateral velocity, percent Initial stream velocity, 0.1 Vpi/il*; temperature, 

porcont ambient, 0.1 V pressure, percent ambient, 0J V mp'^luK - - 

413072—57 15 
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Figure 5. — Disturbances produced behind shook by interaction with 
strongly axisymmetric turbulence. Longitudinal intensity, 0.1 per- 
cent; lateral intensity is 3.61 percent just before shock. Root-mean- 
square velocity in percent of initial stream velocity; root-mean-square 
temperature and pressure in percent of ambient. 

from above, approaching the limiting value of zero. The 
precise asymptotic expression is 

0.^ =0.1 (IJ'VhM-I) 1 " (63) 

where the omitted next-higher-order term is 0[(M— 1)* /4 ]. 

Figure 4 applies when isotropic turbulence flows into the 
shock. Figure 5 (prepared from calculations described in 
appendix B) applies when strongly axisymmetric turbulence 
flows into the Bhock; the specifications for the turbulence 
were taken from theoretical calculations of the modifications 
in initially isotropic turbulence that had passed through 
damping screens and a wind-tunnel contraction (ref. 10, four 
screens, K=2> M=1.5). The calculated deviation from 
isotropy is based on idealized conditions and is pjpbably an 
extreme upper limit to what might be encountered in a wind- 
tunnel test section. The longitudinal component of the 
incident turbulence is the same for both figures — namely, 
0.1 percent of free-stream speed — but the lateral component 
is 3.61 percent for figure 5 against 0.1 percent (isotropic) for 
figure 4. Despite the wide disparity in the lateral compo- 
nent, however, comparison of the two figures shows no great 
change in the curves. Evidently, the lateral component of 
the turbulence flowing into the shock has little effect, and the 



(a) Ambient pressure downstream of shock, 1 atmosphere. 

Figure 6. — Noise generated by shock-turbulence interaction (isotropic 

turbulence). 



(b) Stagnation pressure upstream of shook, 1 atmosphoro. 

Figure 6. — Concluded. Noise generated by shock-turbulence inter- 
action (isotropic turbulence). 

intensity of the remainder of the disturbance field behind the 
shock depends almost solely on the longitudinal component, 
regardless of the degree of anisotropy. The shock-induced 
change in the lateral component itself, however, depends on 
the deviation from isotropy, being appreciable for the iso- 
tropic case and quite negligible for tbe extreme axisymmetric 
case. 

Tbe noise generated by the shock-turbulence interaction 
is measured by the curves of root-mean-square pressure. 
This is best indicated by use of an acoustic scale as in figure 
6. Here the noise pressure level is plotted in decibels above 
tbe standard reference base of 0.000204 microbar for several 
levels of initial isotropic turbulence. According to the pre- 
ceding paragraph there would be little difference for strongly 
axisymmetric turbulence of the same longitudinal intensities; 
the difference between figures 4 and 6 corresponds to no 
more than 4 decibels at the Mach numbers (1.6, 3, and co ) 
for which there are comparable data. 

Tbe reference static pressure behind tbe shock is different 
for the two parts of figure 6. In figure 6(a) the ambient 
pressure behind the shock is constant with Mach number 
(1 atm): this situation may be approximated in an exit jet 
of an aircraft in flight. In figure 6(b) the stagnation pres- 
sure ahead of the shock is constant at 1 atmosphere, so that 
the static pressure behind the shock diminishes markedly 
with increasing Mach number; this situation is roughly 
characteristic of many wind-tunnel flows. It is seen that 
even at a longitudinal component of turbulence of 0.01 per- 
cent, the noise level is severe; and at 1 percent the noise 
level exceeds 130 decibels, which is of the order of the 
threshold of pain, over much of the Mach number range. 

These remarks all refer to the asymptotic noise level an 
“infinite” distance behind the shock, since the attenuating 
part of the pressure waves has been neglected (in practice, 
this distance may be taken to be twice the longest significant 
wave length). For an initial Mach number of 1.5, the noise 
level is predicted to be some 17 decibels greater directly 
behind the shock where the attenuation is nil. 

The local pressure level (proportional to the energy density) 
of the noise field in .the region of shock-turbulence interaction 
is one aspect of the noise problem. Lighthill (ref. 3) has 
investigated another aspect, namely, the flux of acoustic 
energy radiated from the interaction region as a result of 
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the convection of any specified volume of turbulence through 
a weak plane shock segment (1<M< 1.3); the turbulence 
need not be homogeneous. The two quantities, energy 
density and flux of energy, are not simply related unless the 
wave pattern is simple, for example, parallel plane waves 
or concentric spherical waves. 

CONCLUDING REMARKS 

The quantitative effects of the interaction of a convected 
homogeneous field of turbulence with an extended plane 
shock have been calculated, including the pressure level of 
the noise generated in the process. The assumed conditions 
are closely approximated in a supersonic wind tunnel or duct 
with a normal shock: the shock, together with its images 
in the walls (if the latter are nearly parallel), behaves sub- 
stantially like an extended plane shock for the purposes of 
the analysis. The approximation is still quite good for plane 
oblique shocks for that portion of the incident turbulence 
whose eddies are small compared with the tunnel diameter 


(spectral wave length <C tunnel diam.), and probably fairly 
good even without this restriction on eddy size. 

The propulsive free jet emitted by a turbojet, ram-jet, or 
rocket engine is turbulent, but the turbulence is far from 
homogeneous. In addition, only local segments of the shock 
structure that may occur aft of the nozzle can be considered 
sensibly plane. The shock-interaction noise generated by 
turbulent eddies smaller than such shock segments can per- 
haps be estimated from the curves presented herein. Esti- 
mates of this sort refer to the sound pressure level within 
the jet and nearby outside; they provide no direct informa- 
tion on the sound pressure level far from the jet as a function 
of distance and direction, or on the total acoustic power 
radiated by the jet. 

Lewis Fright Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June S } 1954 


APPENDIX A 


SYMBOLS 


The following symbols are used in this report: (In appendix 
B some alternate symbols are defined and used in certain 
parts.) 

a function defined in ref. 2 

a* critical speed of sound 

b function defined in ref. 2 

F(k) arbitrary function of k 

svfr\\ « , , i.* 4a 2 sin 2 0+r 2 COS 2 0 

GY 0 ) screen-effect function, G(B)=— 4 — zr 

w ’ w 4 sin 2 0+m cos 2 0 

27(0) contraction-effect function, 

=Jl 1 

' Zx (e sin* 0 +cos 2 0) 2 

J' Jacobian of transformation from 


e7" 


K 

k 

kr 

k 

dk 


Jacobian of transformation from k to k", 
r// _ dk _ cos 2 6” be' 

J ~dP r "m cos 2 e / W 

pressure drop 


screen coefficient, K- 


dynamic pressure 

amplitude of k: 

radial component of k, k T =—ki cot 6 
wave-number vector, k=k iy k 2y k 3 ; also, 
k=ki, k T , 0 in cylindrical coordinates 
volume element in wave-number space, 
dk—dk x dk>]flk z 

x , final stream speed 

contraction parameter, Zx= . . — -y 

v initial stream speed 

contraction parameter, 

. final stream-tube width 

h — •" 


M 

Mi 


initial stream-tube width 
Mach number-upstream of normal shock 
Mach number downstream of normal shock 


m 


N 

P 


V 

Ri& 

Tfi 

s 


u 

U A 

u , v, w 


V 

V T 

V* 

W 

dWp 


ratio of speeds before and after shock, 

(y+1 )M* 

2+(y-1)M* 

number of damping screens 
transfer function for soxmd waves (pressure 
effect), 

2TymH sec 6 sec O' 

~ (Y+l)m— (y— 1) 
pressure perturbation 
mean static pressure 

perturbation velocity correlation tensor (special 
case of <*/?(£)) 
direction cosines 

transfer function for shear waves, tabulated in 
tables 1 (c) to (k) (eq. in ref. 2 ) 

tan g-l)*+(i tan 6)-, 


Vm (7— 1 )] 

0<^<f 

Transfer function for entropy waves (tempera- 
ture effect) T—Te i5 T 
stream velocity downstream of shock 
stream velocity upstream of shock 
nondimensional disturbance velocity compo- 
nents in directions x u £ 3 , respectively; 

components of velocity perturbation 
’ ’ critical speed of sound a* 

cross-stream velocity (sketch (c)) 
disturbance velocity component in radial di- 
rection/#* 

disturbance velocity component in ^direction/#* 
resultant of U and V 

(complex) wave amplitude in two-dimensional 
field 
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«/S© 

[or/?]© 


transfer function, X=Se is * 

' cos e 

position vector, x=x i =x 1 , 2 ^, 

transfer function, Y—Se n »^^r 

sin 9 

(complex) wave amplitude associated with tv 
(complex) wave amplitude associated with a 
(complex) wave amplitude associated with v 9 
may stand for u , v 3 w, p, p, or r 

screen parameter, for i£>l 

l+il 

correlation of a and p at a separation £ 

Fourier transform of <x/3(£), interpreted as 
spectral density of ai 3 (0) 

rw , 

V cos 2 0' 

perturbation velocity spectrum tensor (special 
case of [a£)fc) 

ratio of specific heats (taken as 1.4) 
phase angle of 1J (eq. in ref. 2) 
phase angle of X and Y, tabidated in tables 
1(c) to (k) (eq. in ref. 2) 
phase angle of T, tabulated in tables 1(c) to (k), 

5 r =tan- 1 (— ^ — \o<0<£ 

\cot 9—aJ 2 

contraction parameter, €=Zf/Zi 

shear-wave inclination ahead of contraction, 

0=tan -1 ^tan 0 ) 

sliear-wave inclination ahead of shock (see fig. 3) 
shear-wave inclination behind shock (see fig. 3), 
0'=tan” x (m tan 6) 

sound-wave inclination behind shock (see fig. 3), 

( , . M* tan 0 n 

—tan 1 — — to-; 0 <9<9ct 
1 —Mr ~~ 

o'— cot" 1 /?., e cr <e<l 


critical value of 9 for which W= speed of sound 
wave-number vector, h~k 1} k 2j k z 
screen parameter, fx=l+a t +K 
screen parameter, v—l+a,—a,K 
separation of two points, 
fimction tabulated in table I (defined in ref. 2) 
transfer function for sound waves, IJ=IIe 15 p 
density perturbation 
mean density 
temperature perturbation 
mean temperature 

perturbation velocity spectrum tensor (special 
case of [a£]®) 

common longitude angle of wave normals 
Jc' y k" in polar coordinates 


transfer function, X=? JJ 


transfer function, T=n 
where 


(cob 9'— e“*” T/a fttt,sin flQ 


(sin 9 / +e inT/2 i 9 W cos 9 t 


n=l, 0 <9<6cr 
=O,0„<0<| 

vf' acute angle between oblique shock and upstream 

flow direction 

Subscripts: 

a, p may stand for u, v , w , p y p, or r 

i, j, m , n may stand for 1, 2, or 3 ; used to replace a and p 

when u, v, w are replaced by u u Uz, re- 
spectively 

Superscripts: 

* complex conjugate 

' refracted shear-entropy wave 

n sound wave 

A distinguishing mark 


APPENDIX B 


COMPLETE VELOCITY SPECTRUM TENSORS 


The first and the sum of the second and third diagonal 
terms of the spectrum tensors of the velocity field behind the 
shock are obtained in the text by use of a simplified approach. 
Other terms are occasionally of interest; for example, the 
separate values of the second and third diagonal terms 
are needed for a description of anisotropic turbulence. The 
complete spectrum tensor for each field (turbulence and noise) 
will be derived herein by a more comprehensive procedure. 

Turbulence field. — It will be convenient to replace the 
symbols u } v, w by Ui, Uz y and to replace a and P by i, j, 
which take on the values 1, 2, and 3 instead of u } v, and w . 
With this notation and the use of equations (8), equations 
(10) can be transformed to 

dZ[=XdZ 1 

dZ^=(Y~l)dZ r cos cp+dZi 


dZi—(Y —l)dZ T sin (p+dZ z 

By introduction of the geometric relations (figs. 1 and 2) 

dZr—dZi tanfl" 
tan d=—kJk T 

f (B1) 

cos <p=—Jc 2 jk r 
sin <p——kz/k T ^ 

all three equations may be represented by the single expres- 


dZ’ { =XdZ 1 S li +\_ (F-l)(-^=) dZt+dZ^il-Su) (B2) 
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where 




fl, i=l\ 

\0, ^F ^ 1 / 


Multiplication of the complex conjugate of equation (B2) 
by the corresponding equation with subscript j and by k* 
yields, after averaging, 

^5Z^=^|Z| 2 ]dZ^5 11 5 </ +(l-5„)(l-5 u )X 

[P l k l kj\Y-l\ 2 \^ 3 +J4dZidZ j + 


hhlcKl -T^dZfdZj+kikMl —Y)dZ^dZf\ + 

(1 - 7 ) [dZ^+kidZfdZA + 
S^l-S^Xlkikikiil-Y^J^+kidZ^dZt] (B3) 

Now, if in equation (15) the symbol for the spectral 
tensor is changed from afi to the more conventional symbol 
ij, application to equation (B3) yields 

ki^tj+kikilciil— Y*)$ lj +k l k,k? r (l — ?)$?<]+ 

Su(l — 5ij)X*\kik]k? T (X — T r )$u+^i , iJ+ 

5 u (l-5„)Z[^ 1 i < F r (l-r*)$ 11 +^$ lf ] y (B4) 

The elements of the turbulence spectrum tensor may be 
exhibited in expanded matrix form: 



It) |X|>*„ 

X> [iihtJa-V) * 11 + 

Y) * 1 + 



t* $»] 

t* <t>iz] 




£)■*»+ 







iihfcj (1-V) 4>n+ 

iibt) (1 -V) *u+ 

1 

I 



Jtibif a-y)*u 




k) <hrl-i5tj|y-l|»tii+ 
fcifaJt) (1- V) 4>n+ 
tibJt) (1-Y) 


The matrix is Hermitian; that is, the missing elements are 


the complex conjugates of the respective elements diagonally 
opposite; that is, and so forth. 

It can be shown, by use of the continuity relation £*$#= 
kj$u=0 (summed over repeated index), that after some 
reduction 

(^+W&'=[(|F| a -i)$„ dk 

in agreement with the second of equations (29). 

Noise field. — With use of equations (Bl), the three equa- 
tions (11) may be represented by the single expression 

dzr t =x dZAt-T (^2) dZ t (1 — Sit) (B 5 ) 

where again the subscripts 1, 2, 3 replace u, v, w, and S u =0 
or 1 as before. Starting with this equation, the spectral 
tensor 3?^ may be derived in a straightforward manner by a 
procedure parallel to that leading from equation (B2) to 
(B4). The result is 

*'tj dk”=$ „ |x | 2 5u5u+|T| 2 ^ (l-5„)(l-5 0 )- 

^5 |X||T|[5 U (1 5u)^+5ij(1 — (B6) 


4^-1 11 % 


(The valid range of this equation has been limited to B„ < 

|0|<^(||| > tan 0 CT ) by use of the simplification xT*=x*T= 

|x||T|, which fails outside that range.) 

The expanded form of equation (B6) is 

|x| 2 # -|x||T|*kW* -|X||T|W 

-|X||T|W-? |T| 2 ^ -|T| 2 Ws 

~\x\\T\kMi -|t| 2 W 3 |T| 2 /r^ 

The diagonal terms yield 

& n dh”= |x|%i dk 

(^+$«)dr=|T| 2 |$„di 

since Jcl-\-Ic 2 =k*; these are in agreement with equations (30). 


APPENDIX C 

CALCULATIONS FOR AXISYMMETRIC INITIAL TURBULENCE 


If the turbulence in the settling chamber of a supersonic 
wind tunnel is considered to be isotropic, by the time it 
reaches the working section it will be axis y mm etric, with the 
longitudinal velocity perturbations very much less than the 
lateral perturbations; the change is due to the effects of the 
damping screens and the contraction (refs. 9, 13, and 14). 
The shock-interaction behavior for a particular case of ex- 
treme axisymmetry will be calculated herein as a matter of 
interest. 

According to reference 14 (with a slight change in nota- 
tion), if the longitudinal spectral density in the settling cham- 
ber (station A') is written as 


[uu]o=k 2 F(k) cos 2 0 (isotropic turbulence) 

then the longitudinal density in the working section (station 
A) is given by 

[uu[=}?F(k) cos 2 06^ (0)i7(0) (axisymmetric turbulence) 

(Cl) 

where k is the wave number at A', 0 is the associated wave 
inclination, N is the number of damping screens, 0(6 ) de- 
pends on the screen pressure-drop coefficient K , and H(6) 
depends on the parameters k and h defining the wind-tunnel 
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contraction. (See appendix A for the functional forms.) 
In what follows, iV=4, 2, Z, =24.92, and 2*=0.3186. 7 

This set of values corresponds (in theory) to an axisym- 
metric turbulence at station A (just upstream of the shock) 
such that the root-mean-square lateral velocity component 
is 36.1 times the root-mean-square longitudinal component 
(see table I, p. 46, ref. 14). The ratio 36.1:1 is clearly an 
extreme deviation from isotropy. 

The effects of the changed form of [uu] on the integration 
procedure will be illustrated by considering the mean- 
square longitudinal velocity in the turbulence. The relevant 
question is (53), with [uu] being given by equation (Cl). 
From the form of equation (Cl) it will be convenient to carry 
out the integrations in terms of x, rather than k ; the trans- 
formation is 


dk=j rj 2 ^£=712 t?d K d<p cos 0 dQ 


Equation (35) then assumes the form 

«*= X rw^' G*H cos 3 ede ( a K*G(H)dK [ iT d<p 
hQjo COS 2 0 Jo Jo 

The last two integrals appear in the expression for izg, the 
mean-square longitudinal velocity at station A r (the ex- 
pression is of the form of eq. (57)) ; thus, equation (C2) may 
be simplified to 


?£=(¥£) f 'Vl 3 GW cos 3 9 dQ 


Uq Wl/Jo 


cos 2 9 


(PS) 


The variable of integration may be changed from 0 to 6 by 
means of the transformation 


1 cos 2 0 
cos 2 6 


dd 


This results in the alternate form 


==(rn) ("Vr Gm cos3 / ?^ e (C4) 
v* W?/Jo cos 2 9 cos 2 9 ^r € 


On numerical evaluation, the integrand of equation (C3) 
was found to have a sharp peak near the upper limit, and 


that of (C4) a sharp peak at the origin. The peaks wore 
avoided by dividing the range of numerical integration 
among the two equations: (C3) was used over the range 
O<0<0i and (C4) was used over the range 6°<0<! 90°, 
where 0i is the value of 0 corresponding to 9= 5°. 
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TABLE I.— WAVE ANGLES 


(a) A/^1.01 (fn=* 1.01669) (abbreviated table) 


dk 

O', 

deg 

8", 

deg 

n 

7.92542 (Bcr) 

a aw 

—81. 944 

1.6346 

7.96000 

a 08i 

-8L 281 

1.3923 

7.97600 

a loo 

-80. ass 

L 3017 

S. 00000 

a 132 

-8a 767 

1.2363 

8.26000 

S.3S0 

-7a 278 

.9132 

8.43786 

a 577 

-7a 471 

.7870 

8.05028 

9.097 

-7a 662 

.6326 

9.46271 

9. 618 

-75.111 

.4633 

9.97514 

iai3S 

-73.630 

.3824 

10. 48767 

10.G59 

-72.339 

.3232 

11.0 

1L 179 

—71. 045 

.2778 

12.5 

12.702 

-67.445 

.1902 

14 

14.224 

-64.014 

.1383 

15 

16.239 

-61.783 

.1142 

20 

20.307 

-60.990 

.0503 

25 

25.365 

-4a 502 

.0251 

30 

3a 412 

-3a 164 

.0124 

35 

35.447 

-19.916 

.0057 

40 

4a 468 

-a 739 

.0010 

45 

45.474 

.384 

-.0001 

50 

6a 466 

10.460 

-.0011 

55 

55.444 

2a 495 

-.0014 

60 

ea409 

3a 493 

-.0014 

65 

65.361 

40.461 

-.0012 

70 

70.303 

5d 401 

-.0008 

75 

75.235 

0a320 

-.0005 

80 

80.161 

7a 223 

-.0002 

85 

85.082 

saiu 

-.0001 

80 

90 

90 

1 


(b) Af-1.05 (w=1.0839S) (abbreviated table) 


deg 

deg 

deg 

n 

16.323(0cr) 

17.61 

-72.39 

1.5246 

16.400 

17.09 

-70.59 

L 2371 

16.630 

17.94 

-03.59 

.9995 

ia936 

ia27 

-easi 

.8376 

17.240 

18. 69 

-05.35 

.7307 

17.649 

18. 92 

—04.05 

.6499 

18. 162 

19.58 

-6L74 

.5336 

18.775 

20.23 

-59.66 

.4608 

19.388 

20.88 

-57.72 

.3870 

20.001 

2L63 

-55. SS 

.3377 

25 

2a 82 

-42.64 

.1357 

30 

32.04 

-3a 76 

.0611 

35 

37.20 

-19. 49 

.0257 

40 

42.29 

-a 62 

.0076 

45 

47.31 

LOO 

-.0016 

50 

62.26 

12.30 

-.0068 

55 

67.14 

22.44 

-.0071 

60 

6L96 

32.41 

-.0067 

65 

66.72 

42.23 | 

-.0055 

70 

7L44 

6L93 

-.0039 

75 

7a 12 

61.53 

-.0024 

80 

8a76 

71.06 

-.0011 

85 

85.39 

sa 54 

-.0003 

00 

GO 

90 

— 


• These values apply for z<=0 only. For values should be replaced by 0. All other? 


TRANSFER FUNCTIONS 


(o) Af-LlO (w«lJ09O8) 


0, deg 

*\deg 

deg 

S 

n 

T 

5., 

deg 

St, deg 

0 

0 

0 

L 145 

-L6667 

-0008792 

0 


180.00 

5 

5.84 

-2a 75 

1.145 

*1. 6503 

-.008788 

.49 

156.00 

10 

1L 65 

-45.45 

1.14S 

•L 6048 

-.008777 

.90 

12a 76 

15 

17.39 

-57.07 

L 153 

•L5318 

-.003767 

1 . 

U 

97.30 

18 

2a 80 

— 61.89 

1.157 

*L 4779 

-.008770 

1.01 

7a 12 

20 

23.05 

-64.51 

1.161 

•L4392 

— . 00S780 

.70 

41.90 

2a 263 

2a 35 

-64.83 

1.161 

•1.4339 

-.008783 

.62 

3a 37 

2a 536 

2a 65 

-65.14 

1.162 

•L42S8 

-.003786 

.51 

29.73 

2a 804 

2a 95 

-65.45 

1.162 

•L4233 

-.008789 

.37 

21.07 

21. 072(0 er) 

24.25 

-65.75 

1. 163 

L4179 

— . 00S791 

0 


0 


21.100 

24. 28 

—64.45 

1.161 

1.25S2 

-.007808 





21.400 

24.62 

-61.02 

1.155 

.9378 

-.005845 





21.800 

25.06 

-sa4i 

1.152 

.7620 

-.004779 





22.200 

25.60 

-6a 34 

1.149 

.6511 

-.004111 





22.560 

25.90 

-54.69 

L 147 

.6769 

-.003663 





24.018 

27.55 

-4a 95 

L 141 

.3S5C 

-.002511 





25.636 

2a 18 

-44.07 

1.135 

.2778 

-.001860 





27.024 

3a 81 

-3a 69 

1.130 

.2066 

-.001428 





2a 612 

32.42 

-35.33 

L 125 

.1561 

-.001037 





30 

34.02 

-3L36 

L 121 

.1187 

—.000373 





35 

39.30 

-ia76 

1.105 

.0448 

-.000373 





40 

44.45 

-7. OS 

1.090 

.0105 

-.000101 





46 

4a 46 

aas 

1.075 

-.0057 

-.000066 





50 

54.33 

14.68 

1.060 

0125 

-.000176 





65 

sa OS 

24.78 

1.047 

-.0141 

-.000252 





60 

ea 72 

34.66 

1.035 

-.0129 

-.000307 





ai 

6a 26 

44.28 

1.025 

-.0103 

-.000346 





70 

72.71 

5a 63 

LC16 

-.0072 

-.000374 





75 

77.09 

62.91 

1.009 

-.0043 

-.000395 





80 

8L42 

72.01 

1.004 

-.0020 

-.000408 





85 

85.72 

8L02 

L001 

-.0005 

-.000416 





00 

| 90 

90 

— 

0 


! ' 





(d) AI-L26 (W-L42857) 


0, deg 

O', deg 

8", deg - 

S 

n 

T 

5., deg 

St, deg 

0 

0 

0 

L 300 

•L 0667 

-a 04059 

0 


180.00 

6 

7. 12 

-laeo 

1.302 

•1.0163 

-.04054 

1.67 

163.61 

10 

14. 14 

-2a 09 

L307 

•1.5840 

-.04010 

2.93 

14a 17 

15 

2a 95 

-3a 66 

1.317 

•1.4882 

-.04018 

4.04 

125.09 

20 

27.47 

-45.31 

1.335 

•L36S4 

-.03998 

4.47 

09.05 

23 

3L23 

-4a 70 

1.353 

•L 2926 

-.04000 

4.11 

7a 76 

25 

3a 67 

-62.33 

1.371 

•L 2466 

-.04019 

3.22 

sa 40 

26 

34.87 

-5as7 

1.383 

•1. 2220 

-.04039 

2.20 

34.30 

2a 300 

35.22 

-6a 93 

L387 

«L 2172 

-.04018 

1.67 

25.44 

2a 660(0 rr) 

35.65 

-64.35 

1. 392 

L 2106 

-.04059 

0 


0 


2a 742 

36.75 

-51.36 

1.373 

.969 

-.03261 





2a 82$ 

36.85 

-50.07 

L 336 

.884 

-.02980 





2a 914 

36.95 

-49.05 

1.360 

.822 

-.02774 





27 

3a 06 

-4a 16 

1.356 

.773 

-.02813 





27.5 

3a 64 

-44.27 

1.338 

.593 

-.02029 





2$ 

37.22 

-41.28 

L326 

.484 ; 

-.01680 





29 

3a 38 

-3a 36 

1.307 

.3484 

-.01237 





30 

39.62 

-32.16 

1.293 

.2614 

-.00953 





35 

4A 01 

—ia 46 

1. 238 

.0658 

-.00270 





40 

5a 16 

-1.86 

1.196 

-.0007 

-.00004 





45 

55.01 

iai3 

1.157 

-.0261 

-.00157 





50 

59.57 

21.02 

1.123 

-.0337 

-.00252 





66 

63.89 

3L.10 

1.094 

-.0327 

-.00316 





60 

07.99 

40.64 

L 063 

-.0277 

-.00360 





65 

71. 92 

49.48 

L047 

—.0211 

-.00392 





70 

76.71 

oa 02 

1.030 

—.0143 

—.00414 





75 

7a 38 

6d 27 

1.017 

— .0084 

-.00430 





80 

82.96 

74.29 

1.008 

— . 0Q3S 

-.00440 





85 

saso 

82.18 

1.002 

-.0010 

-.00446 





90 

90 

90 


— 

— 

* 


' 



are independent of x. 
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TABLE I. — WAVE ANGLES AND 


(e) Af-1.5 (in— 1.86207) 


6, deg 

0*, deg 

0", deg 

S 

n 

T 

i; deg 

Sr, deg 

0 

0 

0 


•1-6667 

-0.1071 

0 


moo 

6 

9.25 

-8.65 


• L 6346 

-.1063 

3.43 


10 

18.18 

-17. 61 

L476 

•3.5441 

-.1060 

a48 


15 

2a 52 

-25.75 

L494 

• L4085 

1047 

8.79 


20 

34. 13 

-33.23 

1.627 

•1.2464 

-.1030 

iaoi 


25 

4a 97 

-4a 01 

1.694 

•L0S50 

-.1019 

9.34 

85.17 

27 

43.49 

— 42.62 

L649 

•1.0366 

-.1030 

7.60 

6b 42 

2S 

44.71 

-43.74 

1.693 

° L0249 

-.1050 

5.38 

4a 29 

28.300 

45.07 

-44.10 

L710 

■ 1. 0240 

-.1058 

4. 10 

29.96 

28.044(0*,) 

45.49 

-44.51 

1.733 

L 0254 

-.1071 

0 


0 


28.670 

45.52 

-42.58 

L699 

.8849 

-.0925 





28.733 

45.59 

-4a 93 

L674 

.7831 

- 0620 





28.760 

45.61 

-4a 60 

1.668 

.7634 

-.0800 





28.822 

45.70 

-39.30 

L 662 

.6973 

-.0733 





28.911 

45.80 

-38.10 

1.636 

• 63S2 

-.0671 





29.0 

45.91 

-37.06 

L 623 

.5913 

-.0624 





29.5 

4a 49 

-32.62 

1.576 

.429 

-.0160 





30 

47.07 

-29.21 

L 644 

.332 

-.0362 





30.5 

47.61 

-2a 31 

1.620 

.286 

-.0294 





31.0 

48.21 ^ 

-23.09 

1.499 

.214 

-.02a 





31.5 

48.77 

-21.30 

1.481 

.174 

-.0199 





32 

4a32 

-19.06 

1.465 

.141 

-.0104 ! 





33 

5a 41 

-14.96 

1.437 

.0905 

-.0109 





34 

51.47 

—lb 21 

1.412 

.0636 

-.0067 





35 

52.51 

-7.73 

1.390 

.0258 

-.0033 ^ 





36 

53.53 

-4.45 

1.369 

.0042 

-.0005 1 





37 

54.52 

— 1.35 

1.350 

-.0123 ! 

-.0017 





38 

55.50 

L60 

1.332 

-.0254 

-.0037 





39 

5a 45 

4.42 

1.310 

-.0356 

-.COM 





40 

57.3S ! 

7.13 

1.300 

-.0435 

-.0068 ! 





45 

6L76 ' 

19.32 

1.231 

-.0618 

— . 0119 i 





50 j 

65.74 

29.87 

1.170 

-.0613 1 

-.0149 





55 

69.39 

39.25 

L 130 

-.0531 

-.0169 





60 

72.77 

47.78 

1.093 

-.0421 

-.0183 





65 

75.94 

55.67 

1.063 

-.0307 

-.0192 





70 

78.94 

63.06 

1.040 

-.0203 

-.0199 





75 

8L81 

7a 09 

1.022 

-.0117 

-.0204 

, 




SO 

84.59 

7a 86 

bOlO 

-.0053 

-.0207 





85 

87.31 

83.47 

b 002 

-.0013 

-.0209 





90 

90 

90 

— 

— 

— 






(g) M~2.5 (m— 3J33333) 


6, deg 

deg 

0", deg 

S 

n 

T 

K deg 

St, deg 

0 

0 

0 

1.700 

• b 0667 

-03139 

0 


18000 

5 

ia26 

-5.95 

b 708 

• b 6S39 

-.3110 

039 

160 66 

10 

3a 44 

-11.86 

1.734 

• b3774 

-.3065 

ias7 

158.51 

15 

4L77 

-17.69 

1.783 

• L 1220 

- 2933 

20.64 

145.33 

20 

6d 60 

-23.43 


-.8638 

-.2722 

2b 12 

127.29 

2a 832 

51.75 

-24.37 


-.8219 

-.2675 

20.80 

123.33 

2L004 

52.94 

-25.31 

b913 

•.7807 

-.2626 

20 60 

118.87 

22.496 

54.08 

-2a 24 

b941 

• 7404 

-.2573 

20 00 

113.76 

23.328 

55.18 

-27.18 

b970 

-.7014 

-.2519 

19.38 

107.69 

24.160 

sa 23 

-28. 10 

2.021 

•.6657 

-.2472 

18.53 

10016 

24.992 

57.24 

-2a 03 

2.086 

" 6375 

-.2448 

17.42 

8096 

25.824 x 

68.20 

-2a 95 

2.209 

-.6323 

—.2511 

15.00 

73.46 

2d 656(0 cr) 

59. 14 

-3a 86 

2.662 

.7638 

-.3139 

0 


0 


2a 675 

59. 16 

—28.88 

2.496 

.6967 

— . 2448 





2a 706 

&a 19 

-27.57 

2.413 

.5120 

-.2108 





2a 756 

69.25 

- 2 a 18 

2.337 

.4366 

-.1801 





2a 856 

6a 36 

-2b 10 

2244 

.3470 

-.1437 





2a 956 

59.46 

-22.61 

2.184 

.2898 

-.1206 





27.056 

5a 57 

-2b 27 

2.138 

.2469 

-.1032 





27.492 

00.04 

-ia70 

2.008 

.1337 

-.0568 





27.792 

ea 35 

-14.29 

1.949 

.0866 

-.0373 





28.328 

60.90 

-ias2 

b 870 

.0290 

-.0128 





28.528 

01. 11 

-a 29 

1.846 

.0129 

-.0057 





28.750 

6L33 

-7.94 

b 822 

-.0026 

-.0012 





2a 164 

6b 74 

-5.63 

b781 

-.0286 

-.0121 





29.264 

6b 84 

-5.10 

b 772 

-.0316 

-.0144 





30 

62.54 

— b40 

b715 

-.0614 

-.0290 





35 

easi 

ia66 

b476 

-.1288 

-.0752 





40 


29.31 

b 340 

-.1276 

-.0932 





45 


3a 24 

b246 

-.1109 

-.1028 





50 

75.87 

47.46 

bl77 

-.0906 

* -.1037 





55 

78.14 

5b 51 

b 126 

-.0705 

-.1126 





60 

sai7 

6074 

b087 

-.0522 

-.1153 





65 

82.04 

6037 

b 057 

-.0364 

-.1172 





70 

83.77 

7b 55 

b 036 

-.0234 

-.1186 





75 

85.40 

7042 

b019 

-.0132 

—.1195 





SO 

86.97 

8b 06 

b003 

-.0069 

-.1201 





85 

88.60 

85.57 

b002 

-.0015 

-.1205 





90 

90 

90 






' 





TRANSFER FUNCTIONS— Continued 


(0 Af—2.0 (m- 2.66607) 



(h) Af-3.0 (m-3.85714) 


0, deg 

0', deg 


S 

n 

T 

S t , deg 

St, deg 

0 

0 

0 

b 741 

•1.0667 

-0.3754 

0 



5 

18.65 

-a 62 

1.751 

•b5610 

-.3725 

11.37 

1©.97 

10 

3b 22 

-lb 22 

1.784 

• b 3125 

-.3631 

ia4s 

159. 10 

15 

45.91 

-10 78 

1.845 

• b 0270 

-.3444 

23.25 

146.05 

20 

6b 64 

-22.27 

b 953 

-.7497 

-.3098 

22.74 


23 

58.58 

-25.62 

2.072 

•.5894 

2707 

20© 


24 

59.79 

-20© 

2.144 

•.5423 

-.2658 

18.56 

97.65 

25 


-27.68 

2.307 

•.5271 

-.2690 

1046 

77.55 

25.3 


-28.00 

2.437 

-.5540 

-.2871 

10 13 

0b 66 

25.644(0*-) 

6b 63 

-28.37 

2.940 

.7137 

-.3754 

0 


C 


25.650 

61.04 

-27.22 

2.785 

.5922 

-.3114 




o 

25.670 

6b 66 

-25. 91 

2.650 

.4872 

-.2504 





. 25.700 

6b© 

-2b 74 

2. 663 

.4118 

-.2170 





25.744 

6b 74 

-23.48 

2.407 

.3456 

-.1825 





25.844 

6b 84 

-2b 41 

2.351 

.2586 

-.1371 





26.944 

6b 95 

-1078 

2.270 

.2040 

-.1086 





20044 


-18.44 

2.221 

. 1649 

-.0882 





20 250 


-10 04 

2.137 

.1072 

-.0579 





205 


-13.61 

2.004 

.0596 

-.0325 





27.0 


-059 

b 961 

-.0019 

-.©11 





27.6 


-5.00 

b875 

—.0465 

-.0266 





28 

wmm 

-3.28 

b 829 

-.0676 

-.0393 





28.5 

6b 48 

-.50 

b780 

-.0874 

-.0520 





205 

M&m 

b 14 

b70l 

1142 

-.07© 





30 

65.82 

027 

b 668 

-.1232 

-.0782 





35 

©.68 

22.74 

b 440 

-. 1489 

-.1180 





40 

72.83 

3b 42 

b314 

-. 1354 

-. 1349 





46 

75.40 

43.68 

b 225 

-. 1130 

-. 1442 





50 

77.73 

fib 15 

bl60 

-.©10 

-. 16© 





56 

7071 

57.63 

1.112 

-.07© 

-.1539 



■ 


60 

8b 49 

63.34 

1.077 

-.0514 

1566 



i 


65 

83.11 

68.48 

1.050 

-.0357 

— . 1584 





70 

8b 61 

73.21 

1.030 

-.0228 

-. 1593 





75 

wmm 

77.04 

b016 

! -.0128 

-.1008 





80 

87.38 

8b 87 

1.007 

-.0057 

-. 1614 





85 

88.70 

86.97 

1.002 

; -.©i4 

-.1618 





90 

90 

90 






H 



° These values apply for only. For x- », values should be replaced by 0. AH other values are Independent of x. 
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TABLE I.— WAVE ANGLES AND TRANSFER FUNCTIONS— Continued 

(I) fmt=»4J57i43) 


(k) Af- co (m- 6.00000) 


$, deg 

tf'deg 

0", deg 

S 

n 

T 

S„ deg 

St, deg 

0 

0 

0 

1.833 

•L0067 

-a 6832 

0 


isa 00 

6 

27.70 

-5.00 

L 853 

-1.4603 

-.6754 

18.07 

17a 58 

10 

40.61 

-iaoo 

1.918 

•1. 0602 

-.5484 

28.05 

ioaio 

15 

68.12 

-15.0 

2.046 

-.7105 

-.4873 

28.51 

14a 67 

20 

65.40 

-2a o 

2.280 

-.3637 

-.3253 

19.33 

moo 

2a 552 

66.03 

-2a 55 

a 317 

-.3140 

-.2889 

17.12 

119.67 

21. 104 

00.64 

-2L10 

a 357 

-.2557 

— . 2420 

14.29 

114.61 

21.050 

07.23 

-2LCC 

a 401 

-.1810 

-. 1761 

10 35 

107.67 

22. 208(0 «•) 

07.79 

-22. 21 

a 449 

.0000 

.0000 

0 


c 


23. 182 

08.73 

-4.93 

1.971 

-.1588 

-.1066 





24.150 

69. 61 

2.45 

1.822 

-.1875 

— . 2004 





25.130 

70.44 

aoo 

L 721 

-.1993 

-. 2301 





26.104 

71.21 

12.78 

L 043 

-.2036 

-.2404 





27.078 

71.94 

lass 

1. 579 

-.2039 

-.2586 





28.052 

72.03 

2a 47 

L 525 

2018 

-.2680 





2a 026 

73.28 j 

23.72 

1.479 

-. 1981 

-.2766 





30 

73.00 1 

2a 69 

L 438 

-. 1936 

-.3819 





35 

7a 61 

3a 83 

L 287 

-. 1638 

-.3021 





40 

7a 77 

47.74 

1.192 

-.1334 

-.3127 





45 

sa 54 

54.76 

L130 

-. 1060 

-.3190 





50 

82.04 

Ca54 | 

1.087 

-.0822 

-.3232 





55 

83.34 

65.48 

1.068 

-.0619 

-.3280 





60 

84.50 

69.82 

L 037 

-.0449 

-. 3281 





65 

85.56 

7a 73 

L 023 

-.0308 

-.3295 





70 

80.63 

77.31 

1.013 

-.0196 

-.3300 





75 

87.44 

80.66 

L 007 

-.0109 

-.3313 | 





SO 

8a32 

8a 86 

1.003 

-.0048 

3318 





85 

89.16 

8a 96 

L001 

-.0012 

-.3321 





90 

00 

90 

— 








a These values apply for j«= 0 only. For x=<», values 3honld be replaced by 0. All other values are independent of x. 


0, deg 

& deg 

0”, deg 

S 

n 

T 

S„ deg 

St, deg 

0 

0 

0 

1.781 

•L 6667 

-04615 

0 


18000 

5 

21.80 

-5.33 

L 795 

*1.5269 

-.4471 

10 95 

170 26 

10 

3a 87 

-10.65 

1.837 

•1. 2248 

-.4327 

20 90 

169.60 

16 

6a 77 

-15.95 

L 917 

*.9087 

-.4029 

25.94 

140 55 

20 

6a 99 

-2L22 

0059 

*.6107 

-.3418 

2036 

127.01 

3L091 

0a 44 

-2a 30 

0100 

* 5444 

-.3204 

21.88 

120. 67 

2i 182 

0L.79 

-2060 

0104 

-.4760 

-.2945 

19.95 

11015 

2a 273 

63.04 

-24.64 

0251 

* 4090 

-.2660 

17.48 

99.12 

23.5 

63.29 

-24.88 

0280 

• 3978 

-.2614 

10 91 

95.02 

24.0 

63.83 

-25.40 

0411 

- 3972 

-.2670 

15.77 

80.18 

24. 304 (^cr) 

04.22 

-25.78 

0283 

.6600 

-.4515 

0 


0 


24.370 

04.22 

-24.47 

0986 

.4707 

-.3281 





24.400 

04.28 

-22.77 

0743 

.3364 

-.2303 





24.404 

04.32 

-20 70 

0553 

.2250 

-. 1545 





24.600 

04.30 

-19.84 

0492 

.1900 

-.1300 





24.604 

04.42 

-1055 

0413 

.1451 

-.1001 





24.064 

64.63 

-1087 

0327 

.0970 

-.0676 





24.760 

04.62 

-15.65 

0273 

.0685 

-.0470 





24.704 

04.63 

-15.45 

0204 

.0641 

-.0446 





25.000 

04.87 

-1067 

0161 

.0111 

-.0078 





25.250 I 

05.12 

-10.20 

0084 

-.0256 

-.0182 





25.773 ! 

65,63 

-5.95 

1.971 

-.0739 

-.0539 





27. 182 

60 93 

064 

1.790 

-.1330 

-.1038 





28.691 

00 13 

9.22 

1.676 

-.1573 

-.1303 





30 

69.25 

14.70 

1.591 

-.1070 

-.1474 





35 

7a 05 

29.37 

1.394 

-.1010 

-.1792 





40 

75.39 

39.94 

L274 

-. 1384 

-.1940 





45 

77.66 

4024 

1.192 

-.1129 

-.2025 





60 

79.60 : 

65.07 

1.136 

-.0890 

-.2078 





55 

SL29 

60.92 

1.093 

-.0678 

-.2116 





GO 

saso 

0000 

1.063 

-.0495 

-.2140 





65 

84.18 

70 68 

1.041 

-.0342 

-.2168 





70 

85.45 

74.94 

1.024 

-.0218 

-.2171 





75 

8065 

7092 

L013 

-.0122 

-.2180 





SO 

87.79 

8071 

1.000 

-.0054 

-.2180 ; 





85 

8090 

8038 

1.002 

-.0013 

-.2189 | 





90 

90 

90 

— 









0) M-QS) (m-5JJ6S29) 



413072 — 37 10 






